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Recently, consensus-type problems have been formulated in the quantum domain. Obtaining 
average quantum consensus consists in the dynamical symmetrization of a multipartite quantum 
system while preserving the expectation of a given global observable. In this paper, two improved 
ways of obtaining consensus via dissipative engineering are introduced, which employ on quasi 
local preparation of mixtures of symmetric pure states, and show better performance in terms of 
purity dynamics with respect to existing algorithms. In addition, the first method can be used in 
combination with simple control resources in order to engineer pure Dicke states, while the second 
method guarantees a stronger type of consensus, namely single-measurement consensus. This implies 
that outcomes of local measurements on different subsystems are perfectly correlated when consensus 
is achieved. Both dynamics can be randomized and are suitable for feedback implementation. 


I. INTRODUCTION 

In quantum as in classical systems, symmetry is a 
fundamental concept, and a key tool in investigating 
dynamical systems. In particular, suitable dynamical 
symmetries, or equivalently the existence of preserved 
quantities, prevent controllability for quantum dynamics 
[aili, while they allow for the existence of protected 
sets of states [UlIiailTlISH]- Symmetric quantum states 
and subspaces have a key role in the description of quan¬ 
tum systems obeying Bose-Einstein statistics [14], they 
are related to thermalization |10j . and have a key role in 
quantum information |21j . Prototypical states used to il¬ 
lustrate and exploit intrisecally quantum correlations be¬ 
tween information units, or entanglement between qubits 
in the quantum information jargon, are the maximally 
entangled states named after Greenberg-Horn-Zeilinger 
(GHZ) [T2] and the W states 0: both are invariant with 
respect to permutation of their subsystems. 

In [20] it is shown that a class of classical dynamics ob¬ 
taining asymptotic consensus can be seen as symmetriz¬ 
ing dynamics, leading to permutation-invariant states. 
The same underlying idea led to dynamics for sym- 
metrizations in the quantum realm m- There, it is 
shown how quantum consensus algorithms can be used 
in combination with simple local controls and measure¬ 
ments in order to prepare pure states and estimate the 
size of a network. This type of symmetrizing dynam¬ 
ics and their convergence properties, as well as their 
continuous-time counterpart, have been further studied 
in P^24j . All the symmetrizing dynamics that have 
been proposed so far, however, are based on combinations 
of permutation operators, and hence they share two com¬ 
mon properties: (1) they are unital, i.e. the maximally 
mixed state is preserved; (2) they attain symmetric state 
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consensus, which is effectively consensus on the statistical 
properties of the variables of interest, but there is no al¬ 
gorithm that can attain actual consensus on the output of 
each local measurement. Due to the contraction proper¬ 
ties of the considered maps, unital dynamics entails that 
the purity of the quantum state can not be augmented 
by the consensus-reaching dynamics. With a degrada¬ 
tion of purity, the fragile quantum correlations encoded 
in the state are typically lost, and stronger notions of 
consensus are out of reach. If the consensus-achieving 
dynamics improve purity, they can be instrumental to 
distributed and robust preparation of interesting states, 
like Dicke (including W) or GHZ states US] [2S|, as we 
shall demosntrate in the following. 

In the effort of overcoming these issues, two types of 
symmetrizing dynamics are proposed, which allow for 
asymptotic symmetrization of the state of a multipartite 
quantum system with respect to the permutation group 
acting on the subsystems. Both dynamics are composed 
of quasi-local maps and are robust with respect to ran¬ 
domization, and are hence suitable for distributed and 
unsupervised implementations - as are their classical and 
classically-inspired counterparts. 

The first one attains an asymptotically symmetric 
state while preserving the expectation value of a global 
observable, reaching average symmetric state consensus 
in the language of m- It does so by selecting and prepar¬ 
ing a pure symmetric state (a Dicke state 13 HD) in each 
eigen-subspace of the observable of interest. It is also 
shown how these dynamics are instrumental to prepa¬ 
ration of globally entangled pure states. The proposed 
protocol uses only single system operations and pairwise 
interactions. 

The second dynamics is the first proposal of a quasi¬ 
local protocol that obtains a stronger type of quan¬ 
tum consensus, called single-measurement consensus m- 
This type of consensus, which implies the symmetry of 
the state but for which the latter is not sufficient, is the 
closest in spirit to classical consensus: after single mea- 
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surement consensus is reached, the measurement of a lo¬ 
cal observable quantity on any subsystem will force the 
whole network of systems to “agree” on the result, i.e. 
yielding perfectly correlating results. Notice that, in con¬ 
trast with classical consensus, the consensus value may 
not be determined before it is actually measured. 

Both methods rely on essentially quantum features of 
the system and, while obtaining a symmetric pure state 
that has a specified average of an observable is not viable 
in general, they allow for a final purity that is typically 
better than the one offered by gossip-type algorithms. 
Both are suitable for implementation via discrete-time 
feedback [3], and can be combined with local initializa¬ 
tion procedures in order to actually prepare perfectly pure 
and entangled symmetric states. This is explicitly shown 
for the first proposed algorithm (see Corollary [^. These 
dynamics can be seen as the discrete-time equivalent of 
conditional preparation of entangled states, in the spirit 
of [26]. 

The paper structure is as follows: in Section]^ a brief 
review of the relevant quantum consensus definitions 
and the ideas underlying existing consensus-achieving dy¬ 
namics is provided; Section |III| begins by explaining why 
trying to obtain quantum consensus with just pure states 
is impossible, and continues by providing the form and 
the convergence proofs for two novel symmetrizing dy¬ 
namics. While these can be straightforwardly extended 
to networks of d-dimensional systems, as those considered 
in the introductory Section, the presentation here focuses 
on qubit networks and pairwise interactions. Qubit sys¬ 
tems are easier to visualize and yet relevant for applica¬ 
tions, some of which have been discussed in [19]. Pairwise 
interactions are the minimal that can be allowed for in¬ 
teracting dynamics, and if a more forgiving locality con¬ 
straint is in place, a set of effective pairwise interactions 


will also be allowed under this locality notion. Section IV 


illustrates the behavior of the two dynamics, and com¬ 
pares them with the existing consensus algorithm. 


II. BACKGROUND 
A. Quantum Consensus States 


erator is the orthogonal complement to its kernel. Given 
a state p and an observable X, the expectation of X ac¬ 
cording to p is computed as Tr(pA). 

For any operator X G denote by A®™ the ten¬ 

sor product X ® X ® ... ® X with m factors. Given an 
operator a G 25(H), denote by the local operator: 

Permutations of quantum subsystems are expressed by a 
unitary operator LQ G il(’H), which is uniquely defined 

by 


C^7r(-^l G . ■ . G Xm)ul = A,r(l) G • • • G 

for any operators Xi,... Xm in ®(H), where tt is a per¬ 
mutation of the first m integers. is the set of such tt. 
A state or observable is said to be permutation invariant 
if it commutes with all the subsystem permutations. 

In [in] a number of potential extensions of the idea of 
classical consensus to a quantum network were proposed, 
and their merit discussed in depth. The ones relevant to 
this work are: 

Definition 1 (SSC) A state p G is in Symmet¬ 

ric State Gonsensus (SSC) if, for each unitary permuta¬ 
tion Utt, 


U.pUl = p. 

Definition 2 (crSMC) Given a G 25("H) with spectral 
decomposition a = ^ •^(^)» « state p G 

is in Single cr- Measurement Consensus (aSMC) 

*/■• 


Tr(nf p) = Tr(nf p), (1) 

for all fc, £ G {!,..., m}, and for each j. 

The definition of crSMC is the unique, among those 
proposed, that requires that the outcomes of a measure¬ 
ments on different subsystems be exactly the same for 
each trial. Consider the set of projections as in 

Definition]^ and let us define 


Consider a multipartite system composed of m isomor¬ 
phic subsystems, labeled with indices i = 1,..., m, with 
associated Hilbert space "H™ := ’Hi G • • • G 'Hm — H®"*, 
with dim(Hi) = dim(H) = n and 2 ^ n < oo. This mul¬ 
tipartite system will act as our quantum network. We 
shall use Dirac’s notation: \ijj) denote vectors of %, (i/'l 
denote their dual linear functional. 25(1^) denotes the set 
of linear operators on %, which in our finite-dimensional 
setting are in a one-to-one relationship with complex ma¬ 
trices. States are associated to density operators, namely 
linear, trace-one, positive-semidefinite operators on 
with their set denoted by 2)(H) C 25("H). Observable 
quantities can be associated to Hermitian operators on 
%, denoted by The support of an Hermitean op¬ 


d 

nsMc = 5In®’". 

It has been shown that a state is in crSMC if and only if 
it holds 

Tr(nsMCP) = 1, (2) 

or equivalently 

flsMCpIIsMC = IIsMCP = P- (3) 

Furthermore trSMC for cr with non-degenerate spectrum 
implies SSC, while the converse implications do not hold. 



3 


Lastly, it is impossible for a state to be trSMC with re¬ 
spect to all a G S){T-L)-. this means that crSMC cannot be 
strengthened to a single-measurement equivalent of SSC. 
The proofs of these statements are given in m- 

It is worth remarking how all these definitions could be 
given for classical systems, in the context of consensus for 
random variables or for probability distributions of the 
state values. In some sense the definition of crSMC is the 
closest to the classical case, as it requires perfect agree¬ 
ment on the outcome of a set of random variables for 
each measurement. However, in contrast with the classi¬ 
cal “perfectly correlated” random variables, the quantum 
version allows for the output to be determined only at the 
moment of the measurement. Both SSC and SMC (with 
non-degenerate observable a) imply that the final state 
is permutation invariant 


B. Quantum Average Consensus 

So far, only “static” properties of consensus states have 
been discussed. However, the core of the problem of inter¬ 
est for this paper is the design of discrete-time dynamical 
systems (or, from an information processing perspective, 
algorithms) that drive the system to consensus. In addi¬ 
tion to this, as in the classical consensus problems, the 
final state will be required to preserve or express some 
property dependent on the initial state. 

Unitary dynamics are not enough when ones is inter¬ 
ested in studying or engineering convergence features for 
a quantum system. A more general framework that in¬ 
cludes (Markovian) open-system evolutions is offered by 
quantum channels [HI [21], that is, linear, completely pos¬ 
itive (CP) and trace preserving (TP) maps from density 
operators to density operators £ : ©("H"*) —>■ S)(H™). It 
can be shown that such maps admit an operator sum rep¬ 
resentation (OSR), also known as Kraus decomposition: 

K K 

^{p)=^ AkpA\ with ^ AlAk = I (4) 

k=l k=l 

where K < (dim(H))^. The representation is not unique, 
however the relation between all the possible different 
representations is well known (see m Theorem 8.2]). 
A CPTP map is said unital if £{I) = I. Given a CPTP 
map £, its dual map with respect to the Hilbert-Schmidt 
inner product £^ is defined by: 

Tr[A£{p)]=TT[£^{A)p]. (5) 

This dual map is still linear and completely positive, 
while the fact that £ is trace preserving implies that £^ 
is unital. 

Locality notions for the quantum network can be intro¬ 
duced as it follows [25| . Consider the multipartite system 
introduced in Section |H A[ An operator U is a neigh¬ 
bourhood operator with respect to a set of neighborhoods 


{AO , j = 1,2, ...,M}, if there exists j G (1,2, ...,M} 
such that: 

V = Vjq-. ® Ijj- ( 6 ) 

where Vjq-^ accounts for the nontrivial action on K.jq'. and 
Ijj) = 0/c^A/; 

Everything is in place to specify the consensus prob¬ 
lems of interest. In addition to asymptotically achieving 
a state with the prescribed “informational” symmetry, 
this state must maintain locally some information on the 
global initial state. Let d{pa,C) = infpgc ||pa ~p||, where 
C C ©("H) and || • || is any p-norm on Given a 

sequence of QL channels {£t(-)}t^o, define £t{po) = Pt = 
£t o £t_i o • • • o £i{pq), and Cssc to be the set of states 
in SSG consensus, and similarly for crSMC. Let S be an 
observable in «B('H®'"). 

Definition 3 (Asymptotic 5'-Average Consensus) 

A sequence of channels {£t(’)}t^o, asymptoti¬ 

cally achieve SSC if 

lim d{£t{po),Cssc) = 0, (7) 

t—¥00 

for all initial states po, and there exists a a G such 

that: 

lim Tr(cr*-^^p(t)) = lim Tr(5'p(t)) = Tr(S'/9o) (8) 

t—yoo t-f-OO 

for all Po and for all £ G {1,... ,m}. 

The same definition holds for crSMC by substituting the 
corresponding state sets in 0. 


C. Gossip-type Dynamics 

Let us recall the structure and core results for the 
gossip-type algorithm introduced in [inj. In a controlled 
quantum network, one can typically engineer unitary 
transformations that implement the “identity” evolution 
and the swapping of two neighboring subsystem states; 
let us denote the latter operator by C/(j j,) for swapping 
subsystems j and k. To develop our analysis, it will 
be convenient to introduce the graph Q associated to the 
multipartite system: its nodes 1,..., m correspond to the 
“physical” subsystems, the edge (j, k) is included if the 
subsystems j and k (have a non-zero probability to) in¬ 
teract. 

Assume edge (j, k) is selected at a certain step t: then 
implement on subsystems j, k the quantum gossip inter¬ 
action: 

p{t -f 1) = £j,k{p{t)) = (1 - a) p(t) + aU(j^k)p{t)Ulj , 

’ (9) 

with a G (0,1). 
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It has been proven that, if the graph associated to pos¬ 
sible interactions is connected, then the quantum gossip 
algorithm ensures global convergence towards: 

= —T Ut,P[)UI € Cssc- (10) 

m\ 

Convergence is deterministic^ when the edges on which a 
gossip interaction occurs at a given time are selected by 
periodically cycling, in any predefined way, through the 
set of available ones; or with probability onel29f. and thus 
in expectation, when the edges on which a gossip inter¬ 
action occurs at a given time are selected randomly from 
a fixed probability distribution {qj^k > 0| k)^E ~ 
1}; Furthermore, S'-average SSC is attained if and only if 
S G can be written, for some a G •Q('H), in the 

form: 

^ m 

( 11 ) 

m ^ 

i=l 

This shows that the mean value of a (global) observable 
can be asymptotically retrieved from the state of any 
single subsystem by employing a quantum gossip-type 
algorithms. 

III. NEW RESULTS 
A. Preliminary observations 


where 0 < A < 1. The only case in which the output can 
be pure for any A is when £[ipk){4’k\) = ^(IV’i)('*Ail) • this 
would satisfy the expectation preservation if and only if 
I'i/'i) and \ip 2 ) had the same expectation, which is clearly 
not possible for j ^ k. Inspired by this observation, we 
will construct dynamics that select a pure representative 
in each of the eigenspaces of the global observable S, yet 
losing some of the correlations between subspaces. 

In the rest of the paper, in order to limit the notational 
complexity and provide some matrix representations that 
are easier to visualize, the focus is on: (1) networks of 
two-level systems or qubits, Tii = span{|0), |1)}. We use 
the common shorthand notation for the standard basis of 
a composite Hilbert space H™: |01...) = |0) (8) |1) 0 ...; 
(2) pair-wise interactions, i.e. the neighborhoods are a 
subset of the set of pairs {{j. A:}, I > j,k < m,j ^ k}. 
Both assumptions can be relaxed, with some straightfor¬ 
ward adaptation in the algorithms. 

B. Symmetrizing dynamics for symmetric state 
consensus 

The preservation of the asymptotic expectations is 
clearly guaranteed if each eigenspace of S is an invari¬ 
ant subspace for all maps that can enter the evolution, 
at any time. More precisely, if S is written in spectral 
decomposition as S' = with H^ the orthog¬ 

onal projectors onto the (degenerate) eigenspaces of S, 
call them "Hfe, invariance of all the Hk is equivalent to 
require 5^(0^) = H^ [5], so that: 


The dynamics associated to the maps are always 
unital, that is the identity is always a fixed point for the 
dynamics. Given the (trace-norm) contraction character 
of CPTP evolutions [1], even if the state asymptotically 
converges to a different state, the evolution is bound to 
get closer to the identity, i.e. the completely mixed state. 

This type of symmetrization, albeit interesting and sat¬ 
isfying the requisites for asymptotic average SSC, is pu¬ 
rity decreasing and in some sense makes the system state 
“more classical”. The natural question is then: can quan¬ 
tum consensus be achieved with a pure state? 

This, in general, is not possible. The desired dynamics 
should select a pure state for any output, while preserving 
the expectation value of the desired global observable: it 
would entail a linear map £ such that: 

Tt{S£{p)) = Tr(5p), (12) 

and Tr£’(p)^) = 1. Write S in spectral decomposition as 
S = J2k^k^ki with lik the orthogonal projectors onto 
the (degenerate) eigenspaces Vk of S. For any pure state 
IV'fc) G Vfe, the output £(\'ipk){'4’k\) = \4>k){4>k\ should re¬ 
main pure. However, any linear map £ would satisfy 

£{M'^k){i^k\ + (1 - A)|^/>j)(V'jj) 

= \k£{Hk){il^k\) + (1 - A)£’(|?/'i)(V'j|) 

= Afe|^fc)((()fc| -f (1 - 


TT{£t{S)p) = Tr(^ \k£t{I^k)p) = Tr(5p). 

k 


On the other hand, as long as convergence to a consen¬ 
sus state is guaranteed, dynamics inside the subspace 
is arbitrary. Let us build one that, while satisfying 
eigenspace invariance, prepares a pure representative in 
each S eigenspace. For networks of two-level systems, up 
to a local change of basis and a rescaling, any S as in 
(111 can be rewritten as 


S = ml+ (13) 

i=l 


This allows to identify the eigenspaces of S with respect 
to this basis as the “excitation” subspaces: 


m 

Tik = span{|n,^2,... € {0,1 },= k} (14) 


generated by linear combination of vectors belonging to 
the standard basis associated to exactly fc I’s. As typical 
in the physics literature, call the number of I’s in such 
vectors, and thus of the generated subspaces, the excita¬ 
tion number. Each of these subspaces contains only one 
pure symmetric state, which is referred in the physics lit¬ 
erature as a Dicke state BE]. Each state is labeled by 
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the number of systems N and the number of excitations 
k, and is associated to a vector: 

\{m,k)) = ]_^0^._0) + ... + I 

y \ k) k m—k m—k k 

(15) 

and each term of the sum corresponds to a unique ar¬ 
rangement of the k ones and N — k zeros. To our aim, 
the key property of Dicke states, beside being pure and 
symmetric, is that they admit a very simple Schmidt de¬ 
composition. Consider a partition of the network into 
two groups of mA and uib subsystems, respectively. It 
is then easy to show that the Schmidt decomposition of 
a Dicke state |(m, k)) is 


\{m,k)) = 


\/(T) 


X! k-kA,kB\{'mA,kA))<S>\{mB,kB)), 

kA+kB=k 

_( 16 ) 

where the weight on each term is HkA,kB = 

In order to explicitly construct the dynamics that pre¬ 
pare these states in each eigenspace, choose a neighbour¬ 
hood Afj, containing 2 qubits, and consider the neigh¬ 
bourhood operator: 




E- 

ieAf, 


(i) 


(17) 


As we did for the global S' in (14), we can construct the 
eigenspaces Sj\f^ as the subspaces spanned by pure states 
with a given number of excitation: 


= span{|ii,z2,.. .,imj),ie G {0,1},^ v = k,} 


e=i 


(18) 


which, in our two-qubit neighborhood, are just three: 


=span{|00)}. 


= span{(|01) + |10))/y2, (|01) - |10))/y2}. 


=span{|l,l)}. (19) 


Joining these sets of vectors (in order), a basis for 
is obtained. With respect to this 
basis, construct a neighbourhood CPTP map £j\f. with 
operators: 



-1 

O 

0 

0 

O 


O 

0 

1 

O 

M,a = 

0 

0 

O 

0 


1- 

O 

0 

0 

O 



1 

0 

0 

0 

Mj 9 — 

0 

1 

0 

0 


0 

0 

0 

0 


_ 0 

0 

0 

I 


( 20 ) 


( 21 ) 


From these neighborhood operators, it is easy to con¬ 
struct a map on the whole network, which is going to be 
the dynamics associated to the selection of neighborhood 
Afj: 

ip) = ® if) (22) 

= MjA®I p 1 © / + Mj, 2 ‘Slip M]^2 S I. 

Here is the identity map on f&i'HjQ-.), which is the 
set of operators that are forced to evolve trivially for 
neighborhood maps on Nj. Now notice that each map is 
CPTP, and: 

1. The first operator is nilpotent and acts like a gen¬ 
eralized “raising” operator, which ensures that in¬ 
variant neighbourhood states for the map £jg'. must 
have support on 

H° = span{|(2,^)),^ = 0,1,2}, 
which is the kernel of Mij- 

2. The second operator Mj ^2 is just the orthogonal 
projector onto The latter is thus an invariant 
subspace for £a(, , s-nd state with support on it 
is a fixed point. 

3. Given the ladder structure of Mij, any state for 
the neighborhood that does not have support on it 
is attracted to Hj. 

In terms of typical available gates/control resources 
for pairs of qubits, such a map can be implemented in a 
“digital quantum simulator” [3], with a two-body uni¬ 
tary (that e.g. maps the ordered standard basis into 
{(|0) + |l))/v^: (|0) - |l))/y2, |00), |11)}), a single qubit 
reset and another two body unitary reverting the effect 
of the first. These maps can also be implemented via 
discrete quantum feedback |1], via a direct projective 
measurement with two outcomes, associated to projec¬ 
tors {Hi — I — Mj 2 ,Il 2 = filj, 2 }i followed by a unitary 
operation when the outcome associated to the first pro¬ 
jector is observed: 


Mjp — HiHi, 


0 

0 0 

o' 

0 

0 1 

0 

0 

I 0 

0 

0 

0 0 

0_ 


where the matrix representation is given with respect to 
the basis constructed above. 

All is in place to prove the following: 

Theorem 1 If the neighborhoods cover the whole net¬ 
works and the associated aravh\30^ is connected, global 
convergence towards SSC is ensured: 

- Deterministically, when all the neighborhoods in any 
order such that all neighborhoods are selected an infinite 
number of times; 
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- In probability and in expectation, when all the neigh¬ 
borhoods are selected randomly from a fixed probability 
distribution {qj = P[A/j] > 0| Qj = !}• of the 

above cases, the system converges to a state p* that has 
support on 


= span{\{m,£)),£ = 0,... ,m}. (23) 


Furthermore, if the initial state has support on the global 
eigenspaces of S with eigenvalue A, the final state is the 
pure state p* = |(to, A))((m, A)|. 


Proof. Let be the span of the Dicke states on the 
whole network. From decomposition (16), it follows that 
for each j, it holds that: 


( 24 ) 


Hence C p|^. 'HP- ® and it is invariant for each 

map as well. It is easy to see that is actually 
= Plj . In fact, if a vector is not in then 

it is also not in for some j. If a state p is invariant for 
all maps , then it has support on the symmetric sub¬ 
space. Since they do not change the excitation number, 
each £jy- leaves the eigenspaces of the corresponding Sjy. 
invariant, as well as those of the global S. 

Consider now a finite sequence of T neighborhood 
maps {£_\f^^^)}t=i, 2 ,...,T in which enough overlapping 
neighborhood are chosen so that they cover the whole 
network of m qubits, and the define the composite map: 


Since each £Mj(t) leaves the eigenspaces of S invariant, 
the same is true for £ff. 

Consider the invariant subspaces 'Hk of S, and define 
14 ™(p) = 1— ((to, k)\p\{m, k)), and the (TO-system) Dicke 
density operators as 


Pra,k = \{'m,k)){{m,k)\. (25) 

Clearly Vf^{pm,k) = 0, and Vf^{p) > 0 if p ^ pm,k, and 
it will be our Lyapunov-type functions: if it is true that 
for all k and p Pm,k with support in Hk there is a 
sequence of maps such that: 

Av^ip) ■■= (C(p)) - vrip) < 0, (26) 


then the only state with support on Hk such that 
= 0 would be Pm.fc- Hence, considering V{p) = 
Sfc ^k^ip)^ main Theorem statement would follow by 
the discrete-time version of the invariance principle (see 
Appendix [a| , since we would prove convergence to H'^, 
and it would be possible to further conclude that Pm,k is 
asymptotically stable in DifHjf). 

The key is thus to prove that there exists a sequence 
such that (26) holds, by induction on the number of 
qubits TO. The inductive hypothesis is: Form qubits and 
all 0 < k < m, there exists £f, such that AVffi{p) < 0 
for any p S D{Hm,k), P 7^ Pm.fc- 


For TO = 2, the network consists in a single pair, so there 
is a single neighborhood Af = {1,2}, and £f = £jy. By 
construction, satisfies the inductive hypothesis. In 
fact, 

H'^ = spanjl 0 , 0 ) 1 , = span{|l,l)} 

only support the invariant states P 2 , 0 tP 2 , 2 , respectively, 
while on H'^^ it prepares the state p 2 .i = (|0,1) -I- 
| 1 , 0 ))(( 0 , 1 | -|- (l, 0|)/2 with a single application, so 
AVfi^ip) = -Fr(p) < 0. 

Next, the inductive step is proved. Choose a partition of 
the TO-l-1 qubits in 2 qubits and the remaining to— 1 ones: 
without loss of generality (up to a relabelling), choose 
Afi = (1, 2} and Afi. Consider a £4-i-i form £_\f^ o 

£P,. Since Pm-i-i k is invariant for all neighborhood maps, 
both V,^+\£Up)) - V,"^+\p) < 0 and V,^+\£Mfip)) - 
< 0 for any initial p [1]. It is next shown that 
at least one of them is strictly negative if p 7 ^ Pm+i,k- If 
it holds already {£tf,{p)) — V^'^^{p) < 0 , then the 

induction step is already proved. Recall that Pm+i,k is 
the projector onto 

span{|(TO-I-1, fc))} = span{|0) 0 |(to. A:)) 

-|-| 1 ) ®\{m,k- 1 ))} 

= span{| 00 ) 0 |(to — 1 , fc)) 

-k(|01)-k|10))(g)|(TO-l,fc-l)) 
-k|ll)) O |(to- l,fc - 2 ))}, 

and 

'H(m-ei,fc) = span{| 0 )}(g)'H(„^,fc)©span{|l)}(g)'H(™,fc_i). 

By the inductive hypothesis and conservation of the 
excitation number k, it follows that £fn{p) guarantees 
< 0 for any state on H{^m+i,k)^ unless: 

supp(p) C (span{|0) (g) |(to. A:))}) © 

(span{|l) © |(to. A: - 1))}). 

Since Pm.fc is invariant for all neighborhood maps, then 
considering f4-i-i(p) = ° ^mip) still makes Vffi de¬ 

crease. 

Consider now the case in which £^ leaves Vffi unvar¬ 
ied and let us consider an additional map £^f„, with the 
neighborhood overlapping, but not being included, with 
the subsystems non-trivially affected by £!f,{p). In terms 
of the Ni,Ni partition, this means: 

SUpp(p) C "Hinv 

= (span{|00) © |(to — 2, k)) + |01) © |(to — 2, A; — 1))}) 

© (span{|10) © |(to — 2, A: — 1)) + |11) © |(to — 2, A: — 2))}) 

The action of £jg'^ maps the subspace orthogonal to the 
target in FAinv, namely the span of: 

14^) = |00)©|(to-I,A:)) 

+(|01)-|10))©|(TO-l,fc-l)) 

— |I1) © |(to — 1, A: — 2)) 
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to 


span {100) 0 |(m — 1, k)) 

+ (|01) + |10))(g)|(m-l,fc-l)) 

-|11) 0 \{m -l,k- 2)).} 

This implies that either p = Pm+i,k or: 

vr^\£MAp)) - vr^\p) = < 0 . 

Notice that is shows it is sufficient to choose the sequence 
of neighbourhoods, and thus the sequence of maps, so 
that they have always overlap with the already selected 
ones and eventually cover the whole network. This is 
also shown by induction on m: for m = 2 it is trivially 
true with no other choices of neighborhoods is pos¬ 
sible. If it is true for m, and the graph associated to 
the neighborhoods is connected, any sequence of neigh¬ 
borhood that satisfies the above properties on some con¬ 
nected sub-graph of m qubits and adding an additional 
two-bit neighborhood that covers the m + 1-th qubit can 
be used. The proof above applies to this case, and the 
Lyapunov function strictly decreases for sequences con¬ 
structed in this way. 

Convergence in probability is then proved by a direct 
application of Lemma in the Appendix, using the se¬ 
quence associated to and convergence in expectation 
trivially follows from convergence in probability. □ 

From Theorem it is easy to prove the following re¬ 
sult. 


Corollary 1 If the following control capabilities are 
available: 

(i) Projective measurements of Oz and unitary control Gx 
for the single qubits; 

(ii) Pairwise interactions as in (22) on a connected 
graph; 

Then any Dicke state is asymptotically preparable. If, 
in addition, projective measurements of S as in (13) are 
available, any Dicke state can be made asymptotically sta¬ 
ble. 


Proof. Consider a target Dicke state, and assume it is in 
the Ceigenspace of S. The strategy to prepare it is the 
following: 

1. First, we want to prepare the correct eigenspace of 
S. To do so, if this additional resource is available, 
first measure S: if the outcome is £, move to the 
next step. Otherwise, if either we cannot measure 
S or the eigenspace is not the correct one, start 
measuring the single qubits. Let I' be the number 
of 1 outcomes, li I — P = q > Q apply Gx to q 
qubits that gave —1 as outcome, if g < 0 apply 
Gz to q qubits that gave 1 as outcome. This will 
prepare a pure factorized, yet not symmetric, state 
in the correct eigenspace. 


2. Next, run the SSC algorithm of Theorem using 
any cyclic repetition of overlapping pairwise inter¬ 
actions that cover the whole network. 

By the previous result, we will have convergence to 
Notice that if S can be measured and 
the initial state is already the correct Dicke state, it 
will remain invariant for the dynamics, as single qubit 
measurements will not be applied. Hence, the above 
protocol effectively asymptotically stabilizes the target. 
On the other hand, if S cannot be measured, the sin¬ 
gle qubit measurements will initially disrupt the desired 
state, which will be re-prepared by the second step, as 
in quasi-local conditional stabilization procedures of [25j . 
□ 


This can be thought as a discrete-time, generalized ver¬ 
sion of the strategies proposed to conditionally stabilize 
W states with continuous semigroup dynamics in [25) . 
Dicke states, with the exception of the ones correspond¬ 
ing to extremal eigenspaces, are entangled and of interest 
for quantum information applications [21] . 


C. Symmetrizing Dynamics for 
Single-Measurement Consensus 


The following dynamics addresses average crSMC. 
Given Theorem 1, it is necessary to design QL dynam¬ 
ics that drive the state on the support of the projector 
IIsmCj while maintaining the expected value of S. In the 
qubit network case, after putting S in the standard form 
(13), the support of Hsmc corresponds to: 


■Hsmc = span{|00.. .0), |11...)}. 


In order to construct the stabilizing neighborhood 

maps, consider the standard basis {|fcj)}^.^o ^ for 'Hjg-., 
where |0) = |00...0), |1) = |00...01), and so on with 
binary ordering until = |11... 1). 

Define = \kj){kj\ and Hjv^^sym = -b 

1 . This is a projector on the neighborhood re¬ 
duction of TLsmCi namely Now construct a 

neighbourhood CPTP map: 

dM„-2 

pNj {p) — IIa/'j .sympHA/"^. ,sym E '^A/j,fc(nfepnfc0?7) 

fc=i 


where: 


T^M,,k{T) = ^ {pkftUkftTUlQ+PkpUkpTUl ^ , 

k=l 

and: 

1. Ukfl satisfies Ukfl\k) = |0) and Ukp satishes 
Uk,i\k) = \dM-- 2 )] 





2. pkfi is equal to the number of zeros in the binary 
representation of k divided by the number of qubits 
in the neighborhood, and pkp = 1 — Pk,o- 

Intuitively, £f^. leaves any state with support on ,smc 
invariant, while on the orthogonal it performs a measure¬ 
ment projecting the state on the standard basis for the 
neighborhood, followed by a “randomizing” map 
This map transfers each basis vector to a mixture of 
|00 ... 0), |11... 1), with weights Pk,o,Pk,i chosen so that 
they preserve the expectation of S. 

The maps on the whole network are obtained by the 
neighborhood maps as: 


(p) = «) /xb (P) ( 28 ) 

By construction, each of these maps is CPTP, and leave 
any state on Hsmc invariant. The following convergence 
result holds: 


Theorem 2 If the neighbourhoods cover the whole net¬ 
works and the associated graph is connected, then the dy¬ 
namics that applies { 28) when the neighbourhood Afj is 
selected ensures global convergence towards aSMC: 

- Deterministically, when all the neighborhoods are se¬ 
lected cyclically; 

- In probability and in expectation, when all the neigh¬ 

borhoods are selected randomly from a fixed probability 
distribution {qj = P[A/}] > 0| % = 1}) independently 

at each time. 


Proof. The density operators with support on 'Hsmc are 
the unique common fixed states for all the maps: if 
state has support on some subspace orthogonal to HsmCi 
there is at least a neighborhood Nj in which the reduced 
state has support outside so it cannot be in¬ 

variant for the corresponding £js/.. Consider 

f4i(p) = 1 — Tr(nsMC,mP), 

where IIsMC.m is the orthogonal projection on 'Hsmc for 
m qubits. It shall be proved by induction on the number 
of subsystems m that there exists a subsequence of T 
maps: 


such that for any p with support outside of "Hsmc: 

i£Vm{p) :=nc(p))-f"(p)<o. 

For m = 2 and £ f = f{i, 2 }: the result is immediate 
given the structure (27): if a state has support outside of 


hirmSMC, it must have support on span{|01), |10)}. As¬ 
sume that q = (01|p|01) -I- (10|p|10) > 0. Then the effect 
of the randomizing feedback map is that of transporting 
the state inside "HsmC) so that: 

AP™(p) := Vm{£^{p)) - V^{p) = -g < 0. 


Now assume that the inductive hypothesis holds for m 
subsystems, with some sequence of neighborhood asso¬ 
ciated to a map If another subsystem is added, the 
effect of £!f^ on this m -I- 1 multipartite system is evalu¬ 
ated. By the inductive hypothesis: 

AK„(p) := Kn(C(p)) - Kn(p) < g < 0 . 

This means that the trace of the output state onto 
'HsMC.m has increased. The subspace 'HsMC.m 0 Hm+i 
has thus also gain probability. This subspace is spanned 
by |00 ... 0) C) lx), |11... 1) ® \y), with x,y G {0,1}. If 
the map £{m,m+i} is applied, anything in this subspace 
is mapped onto HsMC.m-i-i- Thus it also holds that 

AV^+l(p) := Vra+l{£m+li.py) ~ ^ra+l{p) < q <{). 

In this way it has been proved that a sequence of maps 
that contracts Vm for any m does exists. By consider¬ 
ing cyclic application of the sequence associated to the 
map a straightforward application of the invari¬ 

ance principle (in the discrete-time version |18j). asymp¬ 
totic convergence is achieved. 

Convergence in probability is then proved via Lemma[^in 
the Appendix, using the sequence associated to £f„ and 
convergence in expectation directly from convergence in 
probability. □ 

The probability distribution over the neighborhoods 
can actually be made time dependent, as long as all of 
its elements stay bounded away from zero (see e.g. the 
proofs of convergence of [10]). 


IV. A TOY EXAMPLE 


Consider three systems of dimension two, labelled as 
1,2,3, and allow for controlled interactions on neigh¬ 
borhoods {1,2} and {2,3}. The three algorithms for 
symmetrization that are compared, identified by the 
acronyms below, are respectively using: 


(GOS): The gossip-type interactions ([^; 


(DSC): The Dicke-preparing interactions (|2^; 


(SMC): The SMC-preparing interactions (28). 


For all the above, alternating actions on neighborhoods 
Ml = {1,2}, A /2 = {2,3}, are considered. The presented 
simulations illustrate the different behavior of the dy¬ 
namics proposed in the preceding sections. An initial 
state is generated by randomly selecting its spectrum 
and then randomly choosing a unitary orthonormal basis 
in which it is diagonal, and the evolution of the purity, 
along with some relevant state populations, are plotted 
as in Figures |T||3| 

The main observations are summarized in the follow¬ 
ing. 
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FIG. 1: (Color online) Purity and population dynamics for FIG. 3: (Color online) Purity and population dynamics for 

the Dicke states induced by the gossip-type algorithm GOS. the |000), |111) states induced by the SMC algorithms. 



Iterations 

FIG. 2: (Color online) Purity and population dynamics for 
the Dicke states induced by the SSC algorithm 


• Given the theoretical analysis of the previous sec¬ 
tions, the final state for the SSC algorithm has sup¬ 
port only on the span of the four symmetric Dicke 
vectors, and thus the evolution of the four respec¬ 
tive populations is depicted in the figure. Similarly 
for the SMC algorithm and the |000), |111) vectors. 


Figure]^ shows how the average SSC-attaining algo¬ 
rithm associated to maps of the form (22) saturates 
the population of the Dicke states, overall improv¬ 
ing the purity of the state. It is worth noting that 
the increase is not monotone. The natural compar¬ 
ison is with the gossip-type algorithm (|^, which 


also attains average SSC: as illustrated in Figure 
by projecting the state on the group commutant 
muD] with unital operations, the purity decreases 
while the populations of the Dicke states are invari¬ 
ant. While the above behavior is typical for random 
initial states, if the initial states is e.g. pure its pu¬ 
rity will not typically be preserved even for the new 
SSC algorithm, unless it is a state that belongs to 
one of the eigensubspaces of S. 

• The SMC algorithm symmetrizes the state by 
“pushing” it in the subspace generated by 
1000), I 111). The increase of the population of these 
subspaces is depicted in Figure along with the 
associated improvement in purity. A direct com¬ 
parison with the other algorithms shows that, for 
the given initial state, this method leads to the 
highest final purity, with a final value of 0.51 for 
SMC, versus 0.36 for the SSC and 0.17 for GOS. 
This is to be expected: in fact, the SSC dynam¬ 
ics leave all the eigenspace of the operator to be 
averaged invariant, while the SMC leaves invariant 
only the ones corresponding to maximum and min¬ 
imal eigenvalue. Preparing only these two while 
preserving the global expectation allows not only 
for a better purity, but for a perfect correlation of 
local measurement outcomes m- 


V. CONCLUSIONS 

Two classes of QL dynamics that attain asymptotic 
average consensus for networks of quantum systems have 
been introduced. The first improves the existing gossip- 
type dynamics, as in general it attains a purer out- 


25 
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put state while still guaranteeing SSC consensus. This 
is achieved by selecting pure representatives for each 
eigenspace of the global observable whose average has to 
be preserved. While the new algorithm cannot be simply 
obtained by randomized permutations of subsystems, as 
the gossip ones do, it is suitable for a simple discrete-time 
feedback implementation, in the spirit of [3] . The second 
dynamics are the first proposal of dynamics attaining the 
stronger crSMC consensus, thus not only improve the pu¬ 
rity of the output with respect to the gossip-type dynam¬ 
ics, but actually attain a type of consensus that is more 
similar in spirit to classical consensus: if the local observ¬ 
able a is observed for one subsystem and gives a certain 
outcome, it is guaranteed that all other subsystems will 
indeed return the same outcome in subsequent measure¬ 
ments. This second method is also suitable for feedback 
implementation. Both dynamics ensure convergence in 
probability when i.i.d. randomization of the application 
of maps is used, being thus robust with respect to the 
ordering and suitable for unsupervised implementation. 
The i.i.d. requirements can be easily relaxed, following 
e.g. HO]- Remarkably, the proposed SSC-attaining dy¬ 
namics can be used, in combination with few additional 
resources, to asymptotically prepare and stabilize Dicke 
states in multipartite quantum systems. This can be seen 


as a discrete-time version of the conditional stabilization 
protocols proposed in [26] and the results of [T3| for con¬ 
tinuous time quantum semigroups. Similarly, a strat¬ 
egy can be devised to utilize the SMC-attaining dynam¬ 
ics to prepare GHZ-type states with proper initialization 
of the network. As a last remark, the presented con¬ 
vergence analysis is one of the first stability results for 
time-inhomogeneous sequences of CPTP maps using Lya¬ 
punov techniques, and could be of inspiration for more 
general convergence analysis for this important class of 
dynamics. 
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i.e. V{x) is non-increasing along forward trajectories of 
the plant dynamics. Then any bounded trajectory con¬ 
verges to the largest invariant subset W contained in the 
locus E = {a:|AU(x) = 0}. 


Appendix B: Randomization Lemma 

The following Lemma allows us to conclude the con¬ 
vergence in probability of the proposed methods, as soon 
as they admit a finite sequence that ensures a strict con¬ 
traction of the a Lyapunov function. 

Lemma 1 (Convergence in probability) Consider a 
finite number of CPTP maps {£j}fLi: o,nd a (Lyapunov) 
function V{p), such that V(p) >0 and V{p) = 0 if and 
only if p G S, with S some density operator set. Assume 
furthermore: 

1. For each j and state p, V(£j(p)) < V(p); 

2. There exists a finite sequence of maps 

(Bl) 

with j^ G {1,..., to} for all £, such that V{£°{p)) < 
V{p) for all p / S. 

Assume that the maps are selected at random at each time 
t, with independent probability distribution 'Pt[£j\ at each 
time, such that exist some e > 0 for which > £ 

for all t. Then, for any 7 > 0, the probability of having 
V{p(t)) < 7 converges to \ as t converges to -boo. 


Appendix A: Discrete-time Invariance Principle 

The main tool we are going to use in proving conver¬ 
gence of the SSC algorithm is LaSalle’s invariance prin¬ 
ciple, which we recall here in its discrete-time form [15] . 

Theorem 3 (Discrete-time invariance principle) 

Consider a discrete-time system x(t -b 1) = T[x(t)]. 
Suppose V is a function of x G M", bounded below 
and satisfying 


Proof. LaSalle-Krasowsii invariance theorem (in discrete 
time, [T5]) ensures us that by repeating the map £0 in¬ 
finite time would entail convergence to pd, and hence 
V{p) —» 0. Thus, for any fixed 7 , application of the se¬ 
quence (Bl) for a large enough number of times > 0 
would ensure V{£o ^ (p)) < 7 for any p. 

The proof is concluded by noting that the probability 
that a randomly selected sequence oi B ■ N.y ■ K maps 
does not contain contracting subsequence (Bl) is at 
most (1 — . The latter converges to 0 as B, and 

thus t, tends to -boo. □ 


AV{x) = V{r[x]) - V{x) < 0, Vx (Al) 




